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Abstract. In practical work with American put options, it is important to be 
able to know when to exercise the option, and when not to do so. In computer 
simulation based on the standard theory of geometric Brownian motion for 
simulating stock price movements, this problem is fairly easy to handle for 
options with a short lifespan, by analyzing binomial trees. It is considerably 
more challenging to make the decision for American put options with long life- 
span. In order to provide a satisfactory analysis, we look at the corresponding 
free boundary problem, and show that the free boundary - which is the curve 
that separates the two decisions, to exercise or not to - has an asymptotic 
expansion, where the coefficient of the main term is expressed as an integral 
in terms of the free boundary. This raises the perspective that one could use 
numerical simulation to approximate the integral and thus get an effective way 
to make correct decisions for long life options. 

1 Introduction 

1.1 General background 

Initial remarks. The standard model of stock price movement, as proposed 
by Samuelson, is geometric Brownian motion. On the basis of this model, it 
is possible to analyze prices of derivative securities such as options, by use of 
the arbitrage principle, which postulates that a riskless portfolio must, in the 
absence of transaction costs, earn the risk-free interest rate. In this paper we 
shall be concerned with the American put option, which is a contract that allows 
the holder to sell a stock at a fixed price - independent of market movements 
- at any moment during the duration of the contract. It should be pointed 
out that the academic community is not unilaterally in favor of the geometric 
Brownian motion model (see, for instance, 0]). 

The payoff function; geometric Brownian motion. Let f be a time pa- 
rameter, which expresses the time remaining until the option expires. We note 
that t flows backwards with time, so that t decreases as real time passes by. 
Typically, we are interested in an interval < t < T, where t = T corresponds 
to the time when the option is issued and t = is the time of expiration. It is 
convenient to express all money in terms of its equivalent value at the deadline 
t = 0. Let s(t) denote the stock price at time t. If r is the risk-free continu- 
ously compounded interest rate, assumed to be constant, then s(t) exp(— rt) is 
the nominal value (that is, the dollar amount we would see on the screen) of 
the stock at (remaining) time t. For reasons of convenience, we normalize the 
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nominal exercise price of the put option to equal 1. When we think of the stock 
price as a free parameter, we write s in place of s(t). The reward function is 

V r (t, s) — max {0, exp(ri) — s}; (1.1) 

it expresses the payoff earned by exercising the option contract at the point 
(t, s). We consider put options of the American type, which means that the 
holder is at liberty to exercise at any time from purchase at t = T to expiration 
t = 0. Let V r (t, s) denote the correct price of the option at (remaining) time 
t. It is well known from the Black-Scholes analysis that V r (t,s) is obtained 
by optimizing - over all stopping strategies - the expected value of the reward 
function over all stock paths starting at s(t) = s, under the assumption of risk 
neutrality. Risk neutrality means that the expected growth of a risky asset like 
s(t) is postulated to equal that of a riskless one (and since all monetary values 
are discounted to their equivalents at time t — 0, the expected growth is 0). The 
assumption of geometric Brownian motion leads to the infinitesimal equation 

ds(t) = Xs(t)dt + <rs(t)dw(-i), 

where o 2 is the variance per unit of time, A a drift rate (the intrinsic growth rate 
of the stock), and w(— t) is the unit Brownian motion (we write —t to indicate 
that our time parameter flows backwards). The postulate of risk neutrality 
translates into A = 0. By Ito's formula, this then leads to 

2 

d(logs(t)) = y dt + erdw(-i). (1.2) 



1.2 The obstacle problem 

Shift to ordinary Brownian motion. It is more convenient to work with 
ordinary Brownian motion rather than geometric Brownian motion, and so we 
introduce the stochastic process 

. . a y2 , . 

x(f) = — =< logs t , 

V2 a 

which follows a standard Wiener process, as can be seen from l|1.2|) : 

dx(f) = -\/2dw(-t). 
Likewise, we switch from the pair (t, s) to (t, x) as our basic coordinates, where 

X = —= t log s 

V2 CT 

is the free variable corresponding to x(f). In the new coordinate system, the 
reward function l|l.lfl takes the form 

V r , a (t,x) =max{o,e rt -e a2t/2 - CTa;/ ^}. (1.3) 
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Let 

m = 

dt dx 2 

be the heat operator. We say that a function / of the variables (i, x) is caloric 
if 53 / = 0, subcaloric if it is real-valued and ffl / < 0, and supercaloric if it is 
real- valued and S3 / > 0. We need to introduce the upper half-plane (or positive 
time plane) 

R 2 + = {(t,x) : t e R+, x e M}, 

where M is the real line, and R + = {t € 3R : t > 0}. The following statement 
is rather well-known, and may be taken as the formal definition of the envelope 
function V. 

THEOREM 1.1 The envelope function V r ^ is supercaloric in R 2 + . More- 
over, it equals the infimum of all C 2 -smooth supercaloric functions that majorize 
V ri(r (t,x) on R 2 + . 

The scaling properties of the heat operator ffl permit us to reduce the number 
of parameters by 1. We introduce the real parameter g, — 1 < g < 1, and let V e 
denote V r ^, with the parameter settings 

r = l-g 2 , a = V2(l + g): 



V e (t,x) =max{0,e( 1 - e2 ) t - e ( 1 +^ 2t -( 1 +^}. (1.4) 

Likewise, V e denotes the function V r ^ with the same settings. If, however, a 
and r are given, we pick a, g according to 

_ a 2 + 2r _ cr 2 - 2r 

a ^^Mo~' e ~^T2r' 

and recover the function V r ,a from the formula 

V r>(T (t, x) = V e {a 2 t, ax); (1.5) 
by the scaling properties of ffl, we then recover the envelope as well: 

%„{t, x) = V g (a 2 t, ax). (1.6) 

In the sequel, we shall only be concerned with the function V e ; moreover, we 
shall drop the subscript g whenever this does not lead to confusion. 

An affine shift of coordinates. If a function f(t, x) is caloric, then so is the 
transformed function 

e^ t+0x f(t,x + 20t); (1.7) 

here, (3 is a real parameter. The calculation that shows this also reveals that the 
supercaloric functions are preserved under the transformation. Actually, it is 
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possible to give a complete characterization of the transformations of this type 
that preserve the caloric functions. 

Choose /3 = q in the substitution l|1.7ll . and introduce the function 

V' e (t,x) = e^ t+ ^V e (t,x + 2gt). 

Then the least supercaloric majorant to V', denoted by V' is related to V e in 
a straightforward fashion: 

Vfax) = e e2t+ex V e (t,x + 2Qt). 

In other words, we may as well replace V e by V' e in our considerations. The 
function V' is simpler-looking: 

V' g {t, x) = max {0, e t+ex - e^ x } , (t, x)eM 2 + . 

In the sequel, we shall consider only the transformed function V' and write V e 
for it. 

Introduction of a new parameter. We introduce the parameter 9, confined 
to < 9 < +oo, and let V g> $ be the function 

V e , e (t, x) = e t+ ^ - 1(1 + 9 - g + 9g) e*"" - 1(1 -0)(1 + g) e t+x , (t, x) G n +! 
extended to the whole positive time half-plane by 

v e , e (t,x) = o, (t,x)eR 2 + \n + . 

We readily calculate that on Ml , 

ffl V e , (t, x) = -9(1 + g) e* S (x) + (1 - g 2 ) e t+ex l n+ (i, x), (1.8) 

where 1^ denotes the characteristic function of the set E, and <5o is the unit 
Dirac mass at 0. 

From this, it is evident that the role of 9 is to scale the mass distribution 
along the t-axis. The value 9 = 1 corresponds to the put option problem in the 
introduction. We have the corresponding envelope function V e ,g, and Theorem 
II. II generalizes to the new setting. 

1.3 The free boundary 

The continuation region. The region 

V(g,9) = {(t,!)?^ : V g , e (t,x) < V e ,e(t,x)} 

is called the continuation region, while the boundary curve 

T(q, 9) = dV(g, 9) n 
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is the free boundary, or decision boundary. Note that 

R+\n+ c v(g,e). 

We parametrize T(g, 9) by 

£ = < i < +00, 

and note that it is well-known that (f> e ,e{t) is an increasing function, with 
05,6(0) = and 



km </> Bi g(t) = fi( e , 9) = -i- log (l + 



(1.9) 



where the equation is taken to define the constant ji = fj,(g,6). Moreover, at 
least for 9 = 1, it is known that the function a; = 4> g ^e{t) is continuous and 
concave (see, for instance, [H]). In financial terms, in the continuation region, 
we keep the options contract, while at the decision boundary, we exercise it. 
The behavior of the free boundary near t = has been studied extensively (see 
HI): 

<f>e,x{t) ~ (l + f?)^2tlogi i^O (for 0=1). 

Here, we shall focus on the behavior of the free boundary as t — > +00. We shall 
demonstrate that <p g< e (t) has the asymptotic expansion 



where /3j = (3j(g,9) are certain real- valued coefficients. The term 

stands for an L p (M + )-function times i~ Ar_1 e _t , and we are free to choose p 
with 2 < p < +00. We also show how to express (3j in terms of an integral, 
which should permit numerical algorithms to yield good approximations of (3j . 
In particular, 

p —gH f+00 

/?! = —/ <jy e . e (t) e°*-» w 4>' ej g (t) e* dt, 
2 V 7r Jo 

which we may - by integration by parts - also write in the form 
1 r<9/ie-^ (1 + gn)e- e ^ -1 

/3l = ^\T^ + J 



1 

Q Jo 



+00 

' {QH - Q<t>eAt) + l)e e ^.«W-^ - 1 e* dt 



since it is a consequence of our analysis that 

"+00 



<f>' e fl (t)e e e * d* = 
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1.4 Asymptotic estimation of the free boundary 

Estimation of V e fi from above and below. We need to get a grasp of the 
asymptotic behavior of the function 4> et $(t) as t — > +00. Let V e fi be the solution 
to the heat equation with boundary values equal to V g .g on the two half-lines 

{(i,i)6i 2 : i = 0, -oo<x<y) 

and 

{(i,i)£l 2 : < t < +00, x = 

it is unique under mild growth restrictions. This function arises from the stop- 
ping strategy of stopping on the two given lines, and hence we must have 

V e ,e(t, x) < V B>e (t, x), (t, x) e n M , 

where 

n M = {(t,x) G R 2 : 0<t<+oo, -oo<x</j,}, (1.10) 

and // = fj,(g, 9) is as in l|1.9fl . 
Let Voo be the function 

Voofax) = ve t+x , {t,x)eR 2 + , 

where 77 is the constant 

^ = i(i + e)[ e (»-D/'-i + fl]; 

then 

Vg,e(t,x)<Voo(t,x), (t,x)eU„. 
A calculation reveals that 

V 00 (t,x)-V g , e (t,x) = -{l-Q^e^etix-tf + O^ix-n) 3 ), (1.11) 

near the line x = \i. Let V\ denote the function 

Vi(t,x) = Voo(t,x) - v Bt e(t,x), (t,x) e n M , 

and note that it is caloric and that it vanishes along the boundary vertical line 
x = [i. By the reflection principle for caloric functions 5, pp. 115-116], then, 
V\ extends to a caloric function throughout the positive time half-plane, with 

V 1 (t,x) = -V 1 (t,2fj,-x), (t,x)eR 2 + . 

The boundary values of V\ along the a;-axis are given by 

Vl (0, X) = Voo (0, x) - V e>0 (0,x) = Voo (0, x) - V e , g (0,x), X &] — CO, fj]. 



6 



This is not as explicit as desired, so we decompose 

V 1 (t,x) = V 2 {t,x)-V 3 {t,x), 
where V2 and V3 are caloric in the positive time half-plane, with boundary values 
V 2 (0,x) = V sgn(p - x) min {e x , e 2M ~ x }, 

and 

y 3 (0, x) = V{x) = e ex - 1(1 - e )e^ e»- x - 1(1 - 9)(l + g)e* , x e [0, M ], 

with V3(0, x) = for x e M \ [0, 2/z], and V 3 (x) = -V 3 (2fj, - x) for x € [m, 2/z]. 
The first boundary moment of V2 is 

/+00 
(a;- M )V 2 (0,a;)dx = -2r7, 
-00 

and that of V 3 is 

(x - M ) V 3 (0, x) dx = -2 

-00 (? 

+ (1 - Q )eS e+1 ^(e-» - 1 + + (1 - + e )(e" - 1 - m), 

so that 



/+00 
(a; - /u) Vi (0, z) da; = — 2r? + 2 
-OO 



- 1 - (?m 



-(l-e)e( e + 1 )^(e^-l + /i) + <0; (1.12) 

after all, the function V\(Q,x) is positive for x < fi and negative for x > \l. We 
need the following observation. 

LEMMA 1.2 Suppose f is a continuous odd complex-valued function on the 
real line, which decays exponentially rapidly at infinity. Then its caloric exten- 
sion f(t,x) to the positive time half plane, as given by the formula 

1 f +oa (x-0 2 
f(t,x) = -= / e"^/(0d£, 

ZVTTI J -oo 

has the asymptotics 

as t — > +00 and x is kept inside a compact interval of the real line. 
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The proof is left as an exercise to the reader. 
In view of the lemma and we obtain 




so that V e fi{t, x) < V e _g(t, x) holds in a domain of the type 



x-n< -Bit- a/2 e-* + 0(t~ 5/2 e- t ), 



where the positive constant B\ is given by 



Bi 



20F(1 - Q 2 ) 



e 




(1.13) 



which is evaluated in equation l|1.12|l . We have obtained the following statement. 

LEMMA 1.3 The function 4> e ^ describing the decision boundary enjoys the 
following estimate: 

[i - B 1 r 3/2 e- f ' + 0(t- 5/2 e- f ) < <j> et g{t) < fi, < t + oo, 

where the positive constant B% is given by fl.lty . 

Note that this proves the assertion that the vertical line x = fi is an asymp- 
tote for the decision boundary. 

2 The balayage equation 

2.1 Derivation of the balayage equation 
Integration by parts and balayage. We introduce the function 



in the continuation region V(g,9), it solves the overdetermined problem 



U e fi(t,x) = V e ,e(t,x) - V e ,e{t,x); 




o 



on V(g,6), 
on &D(g,6), 
on T{q,6). 



In the complement M. 2 , 



we have 



V e ft{t,x) = V e ,e(t,x) 



(t,x)eR 2 + \V(g,9), 
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which makes U B: g — there. Also, as U Bj g = on &D(g,6), we can extend 
the function U Bt g continuously to all of R 2 by declaring V = throughout 
R 2 \ T>(g,6). Since U e .g vanishes together with its gradient along F(g,9), it 
actually solves 

f EB U e . g = -lv {e ,e) ffl V e , e on R 2 , 
\ f/ e ,e = on R 2 \V(q,6). 

Written out more explicitly using 1)1. 8fl . (|2.1[l assumes the following form: 

fflC/ e , e = [0(l + g)5 o (a;)-(l-g 2 )ln + nD( e ,e)(i,x)] e t+ex , (2.2) 



on R^_ , while 



U g .g{t,x) = 0, (t, i)cl 2 \ T>(g, 9). (2.3) 



By integration by parts, if h is a compactly supported C^-smooth function in 
R 2 , then 

/ SU e ,g(t,x)h(t,x)dtdx = U e , 9 (t,x) m*h(t,x)dtdx, (2.4) 

where 

dt dx 2 

is the adjoint heat operator. We would like to plug the functions 



h(t, x) = h z (t, x) = e 



-tz +XZ 



into (|2.4I) . for z E C, because they are all EB*-caloric. This function, unfortu- 
nately, is never compactly supported; however, it might be possible to approxi- 
mate it by compactly supported functions so that l|2.4|l holds for h = h z in the 
limit. As we proceed in this manner, taking into account the known growth 
properties of V g ^, we find that we should restrict the complex parameter z to 

1 < Re z and 1< Re(z 2 ). (2.5) 

We compute 

(1 + q) 5o(x) e t+ ** e~^ +xz dtdx = ^1±|> , 

z z — 1 

■ ~+ 
while 

(1 - Q 2 ) ln + nv( e .e)(t, x) e t+ex e~*^' dtdx = 

i 2 /-+oo -i 2 



Q- z Jo (g + z)(z 2 - 1) 



9 







It now follows from (|1.8|l and (|2.4[) that 

We shall call this the balayage equation for the free boundary equation x — 
4> e ^e(t). The reason is that the positive mass concentrated on the half- line 

{(t,x) eR 2 + : x = 0} 

should be counterbalanced by a corresponding negative mass spread evenly over 

u+nv( g ,e) = {(t,x) el} : < x < <t> e , e {t)}- 

It is possible to show that, under reasonable restrictions, the balayage equa- 
tion characterizes the free boundary. We shall, however, not pursue this matter 
here. 

2.2 Equivalent formulations of the balayage equation 

Rewriting the balayage equation. In the context of the balayage equation 
Q2.6p . we introduce the complex variable s = z 1 . We use the principal branch 
of the square root to define y/s, so that if s has positive real part, then so does 
yfs. In particular, if z meets (|2.5(l . then y/s = z. It follows that we may rewrite 
H2.6JI in the form 

+ °° e e , e (*)[ e +vq e -[*-l]t dt = 1-Q + 0(Q+^) ^ (2 ?) 

for all s € C with Re s > 1 . 

We introduce the function fg^e-, 

ip e ,e(t) = ft — <j> et e(t), ie[0,+oo[, 

where fi = fj,(g, 9) is as before. The function (p g ^ is decreasing, with 

f e ,e{0) = M 

and 

<p efi (t) = 0(^ 3/2 e"*) as t -> +00, (2.8) 

in view of Lemma 11.31 In terms of this function, the balayage equation (|2.7|l 
assumes the form 



7o 1 - £> s - 1 



(2.9) 



for Res > 1. This is the version we shall use many times in the sequel. By the 
way, an integration by parts manoeuvre applied to l|2.7|l results in the simpler- 
looking identity 
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valid for Res > 1. In terms of the function ip e $, the relationship reads 

k eW|e" Ve ' 8(t) L e+yi J e-^-^dt= — e - ^, (2.10) 

y ' 1 - £1 

where, again, Res > 1. 

2.3 A general scheme for analyzing the balayage equation 

Taylor's formula for the exponential function. The exponential function 
has a Taylor expansion about the origin: 

z 2 z N 
e z = l + z+ — + ...+ — + E N+1 (z), 

where the remainder term En+i{z) can be expressed in the form 

1 r z 7 N+1 r 1 

E N+1 (z) = — J J [1-tfe^dt, 

which represents an entire function, and enjoys the estimate 

\E N+1 (z)\ < JJ_^ m ax{l, e R -}. (2.11) 
As we apply the above formula to the balayage equation l|2.9|l . the result is 



s 

\N „, r+oo 



+ m~ [Q + V5] N C[<t$, ] {s-l)+j' E n+i (-[d+VZ\ <P e ,0 (t)) dt 

'^l ll±!Me ^ 

1 — Q S — 1 

for Res > 1. Here, C denotes the Fourier-Laplace transform, as defined by 

C[f]{z)= / e~ tz f(t)dt, 
Jo 

wherever the integral converges. The intention is to use the identity (|2.12l) to 
sucessively obtain more information regarding the function tp e fi- Indeed, in view 
of Lemma fl .31 we have some input to initiate the iterative process. 

2.4 Asymptotics of the free boundary 

Analytic continuation of the Laplace transform. We recall the definition 
of the function tp, 

4>g,e(t) = M - ^gfiii), < t < +oo, 
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which is a decreasing function with tp g> g (0) = fj, and the asymptotic bound 12. 8|) . 
We use the representation l|2.12[l with z = y 7 ? and N = 1 in the form 

r\,n u,a n 1 e ~ W l- e + fl(z + e) __ M 



z 2 - 1 1 - g 



E 2 ( - [q + z]<p e , e (tj) e- (z2 -^ dt. (2.13) 



It is immediate from the definition of /i — fi(g, 0) that the function 

Z Z — 1 1 — Q Z Z — 1 

extends analytically to C\ {— 1}, and that the singularity at z = —1 is a simple 
pole. Moreover, by Lemma 1 1.31 the left hand side converges to a holomorphic 
function for z with Re[z 2 ] > 0. Furthermore, by H2.11(l . the function 

r+oo 

S 2 (z) = J E2[-[g + z}ip e ,e(t))e-^ -^dt 

expresses an even analytic function in the simply connected domain Re[z 2 ] > — 1, 
which extends continuously to the closed region Re[z 2 ] > -1. It follows that 
the right hand side of l|2.13|l expresses an analytic function in Re[z 2 ] > — 1 with 
the exception of a simple pole at z = —1. 
We introduce the notation 

<J>(z) = £[^, e ](z 2 -l) 

for z £ C with Kez > and Re[z 2 ] > 0, and let <& denote the possible analytic 
continuation of this function beyond its initial domain of definition. Then H2.13J1 
reads 



1 e-w l-g + d(z + g) 



Z 2 -l l-Q Z 2 -l 



E- 



(-[g + z^^t^e-^-^dt, (2.14) 



for all z € C with Re[z 2 ] > —1 and z ^ —1, and $(z) is holomorphic in 
Re[z 2 ] > — 1 with the exception of a simple pole at z = — 1, and continuous up 
to the boundary. 

Estimate of the Laplace transform. We need some size control of &(z), in 
order to invoke inverse Laplace transformation and obtain information about 
<Pq.8- By integration by parts, 



z—J H-^ + ^At))e^-^dt = ^- 1 



z + g 



<ff Bfi {t) E x ( - <p B , e (t)[* + Q\) z~ [z2 ~ 1]t At - ( 2 -15) 
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Combining l|2.13[l and (|2.15l) . we arrive at 

a P -M(z+e)" 



z 2 - 1 



1-f? 

1 



#i ( - + Q]) e- [z2 - 1]t dt. (2.16) 



z 2 - 1 

Let us, as a matter of convenience, introduce also the function 

*(z) = n - [z 2 - 1] $(z), 
which is a Fourier-Laplace transform as well: 

*{z) = C[- l p l efi ]{z 2 -l), 
for z € C with Rez > and Re[z 2 ] > 0. In terms of 12.16|l reads 

*(*) = — — + / V\t) Ex ( - <p e ,g(t)[z + g}) e-I* dt, (2.17) 



for z S C with Re[z 2 ] > —1. By letting z — > along the positive real axis, we 
obtain from l|2.17|l that 

*(0)=/ |¥/ e)fl (t)|e t d*=- / \< f / Bi g(t)\Ei(-<P*At)Q)e t to> 

Jo 1 — Q Jo 

which simplihes to 

y e ,e{t)\e-^^e t dt=- . (2.18) 

1 - Q 

From this and the fact that the Laplace transform of a positive function is 
dominated by its behavior along the real line, we get, by integration by parts, 
that 

|£[<p e ,e](s)| = °(j^J as M ->■+«>, Res > -1. (2.19) 

This leads to an estimate of $(z) in the quarter-plane Rez > 0, Re[z 2 ] > 0. We 
need size control of &(z) in the bigger region Re[z 2 ] > — 1. From 1)2.17)1 and the 
estimate 1)2.11(1 . we obtain 

flp-M(e+Rez) 

|*(z)| < 

1 - Q 

r+OC 

+ + / \^ e (t)\ V Mt) ™*{he~ v ^ t)le+Rcz] } e 2t dt. 
Jo 

In view of the derived integrability properties of ip and its derivative, this gives 
the following growth estimate for <E>: 



Mz)\ = o 



X max{l,e^ Rcz } ) , as \z\ -> +oo, Re[z 2 ] > -1. (2.20) 
\z\ J 
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The Laplace transform and a slit domain. It follows from the derived 
properties of $(z) that the Laplace transform C[ip e ^](s) extends analytically to 
the region Res > — 2 minus the slit ] — 2, — 1] on the real line, and it extends 
continuously to all boundary points with Res = —2, except for s = — 2. At 
s = — 1, C[(pg t g](s) has a square root branch point, and for real x, —2 < x < — 1, 
we can speak of the functions £[ip e> g](x + iO) and £{<p e j](x — iO) as continuous 
limits from above and below. This makes the following expression well-defined: 

AO) = ±.(£[<p Btg ](- x + iO) - C[(f e ,e](-x - iO)), 1 < x < 2. 
Using the identity (|2.16|) . we can write this function as 



A(x) 



irx I 1 — Q 



/XV 'x — 1 



1 < x < 2, (2.21) 



from which we see that it is real-valued on the interval in question, and has a 
square root type singularity at the point x = 1. We now extend the function 
A to the whole real line by setting it equal to off the interval [1,2]. We shall 
need the function ip\, the Laplace transform of A: 



VA (t) = £[A](t) = 



-Hoc- 



A(x)e- tx dx. 



It is real-valued, and, due to the noted property 



A(x) = 0(Vx~l) 
it has the following asymptotics: 

The Laplace transform of ip\ equals 

C[ip g .e](x + iO) - C[(p g .g]{x - iO) 



as t 



£[va](s) 



1 

2tt7 



dx, 



which means that the difference 



C[tp efi ](s) - C[ip A ](s) 

is holomorphic throughout Res > —2, with continuous boundary values, except 
possibly for a logarithmic singularity at s = — 2. It has the decay rate 

|£M(s)| = O (1) , \s\ -> +oo, Res > -2, 
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so that according to (|2.2U|) , we have 

\C[ip — (pA_](s)\ = O ( — ) , as |s| — > +00, Res > —2. 



1*1, 

Using some basic complex analysis, we get the intermediate growth control 

\£[cp-cp A ](s)\ = o( ^-^^ j , as|s|^+oo, -2<Res<-l. 

In particular, by the Plancherel identity, the function 

tH->e< a - e >*(w(f)-p A (f)) (2.22) 
is in L 2 (R + ) for each positive e. Together with 

<P e ,<?(*) ~ <PA(t) = 0(t" 3/2 e-') asi->+oo, (2.23) 

this gives us rather good asymptotic information regarding the behavior of the 
difference function tp Bt g — ip\. 

2.5 The asymptotic formula 

Generalized Taylor series for A. The function A(x) has a convergent series 
expansion 



A(x) = y/x-l^Xo + Xi(x - 1) + X 2 {x - l) 2 + . . . ) , 



1 < x < 2. 



It is well-known in asymptotic analysis (see PJ) that this leads to an asymptotic 
formula for its Fourier-Laplace transform tpA, 

Mt) = E a, r (j + §) / ' ■ \ ' + o(^- 5 / 2 e-*) 

as i — > +00. In view of lj2.22|l and H2.23(l . we get the same type of expansion for 

<p,e (t) = J2 A; r (j + I ) i^" 3 / 2 ^ + O iP (t^e"') , (2.24) 

where by Oj>(l) we mean an expression with bounded norm in L P (R + ), and 
2 < p < +00 is arbitrary. The coefficients Xj may be read off from the identity 
12.21fl : for instance, if we use in addition the identity l|2.18|l . we obtain 

fill ^~M£? 1 /' + 00 



7r(l — Q) 7T 



^, e (t)|(M-^W)e-^« e e*dt. (2.25) 
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By QFZty with N = 1, this leads to 



<P e ,e («) = Pi t-^e-' + O lp (t- 2 e- ( ) , 



for any p, 2 < p < +00, with (3\ = T(3/2) Ao- This is the formula alluded to in 
the introduction, subsection II. 31 

3 Further topics 

Maximal meromorphic extensions. It is of interest to apply the general 
scheme <|2.12ll not just to N = 1, but also N = 2,3,4,.... Unfortunately, 
this approach seems to run into difficulty after a few steps. In any case, it 
is possible to show that C[ip Bt e]{s) has an analytic extension to the half-plane 
Res > —4 minus the slit ] — 4,-1] along the real axis. One possible approach 
to get an explicit expression for C[ip Bt s\ (s) would be to obtain the maximal 
meromorphic continuation to a Riemann surface sheeted over the complex plane. 
If the structure of the Riemann surface would happen to be simple, it might be 
possible to express C\<p Qt g\(s) in terms of the poles of the maximal extension. 
This would then also give an explicit way to compute the free boundary function 
<p e ,e(t) itself. 

Small values of 6. The parameter 9 was introduced mainly for the purpose 
of having the opportunity to choose close to 0, since everything is completely 
understood for 9 = 0. It might then, in a second step, be possible to extend the 
analysis to general values of 9. 

As we differentiate the balayage equation l|2.9|l with respect to the parameter 
9, we have 




d// 



1 



+ {Q + Vs) 



9 dfi] e-^+v^) 



1 



1 



g d9\ s-1 



and after another differentiation, we have 
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For 6 — 0, <fi e ,$(t) = 0, and as we apply the first identity, we realize that we also 
have 



dd 



(*) 



= 0, < t < +oo. 



Inserting 6 = also into the second, we get 



d A <p 



(*) 



1 



1 



d<9 2 

and, as a consequence 

d 2( Pg,e 



e -(s-l)t dt _ 



Res > 1, 



d(9 2 



■(*) 



e=o 2(1 



-1/2 /-+QO 

^ / r- 3 / 2 e- T dr, < t < +oo. 

(1 - Q) 2 Jt 



This gives us some understanding of the behavior of <p B ,9(t) for small 6 and fixed 
q and i. 
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